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We investigate QCD sum rules for vector currents in the rho meson channel in the nuclear 
medium. For increased sensitivity, we subtract out the vacuum contributions. With a saturation 
scheme often considered in the literature, we find these "subtracted" sum rules to be unstable. It 
indicates the importance of other contributions neglected in the saturation scheme. These include 
' more Landau singularities at low energy, additional operators in the medium, and possibly the 

f^*) , in-medium width of the rho meson. 
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^ ; I. INTRODUCTION 

£C) , It was quite some time ago that the original vacuum QCD sum rules M were first extended to a medium at finite 
temperature and density [^.^). Since then much work has been devoted to extracting results from them. Typically 
one tries to predict the medium dependence of the parameters of particles and resonances in the low energy region. 
However, these in-medium sum rules have several new features. These must be included quantitatively in the saturation 
' scheme for reliable predictions. 

One such feature, recognized soon after their initial formulation, is that the in-medium sum rules have contributions 
^vq from new operators over those already present in the vacuum sum rules ■ The velocity four- vector of the medium 
allows one to construct additional Lorentz scalar operators. (In the rest frame of the medium they are simply the 
t-H • time components of Lorentz tensor operators.) At higher dimensions there arises a multitude of such operators. 

The other features relate to the hadron spectral function. A communicating single particle state may acquire a 
• large width (in addition to a possible shift in mass) due to its scattering from the particles in the medium, so that it 
|' may not be possible to treat it in the narrow width approximation as in the vacuum. Also multiparticle states in the 
JL , medium may give rise to several nearby singularities of the Landau type. 

In this work we reexamine the sum rules in the p-meson channel in the nuclear medium at zero temperature. Our 
saturation scheme is similar to one mostly employed in the literature, which includes only some of these features 
In the operator product expansion we include properly all the operators up to dimension 4 and only the four-quark 
. £^ \ operators at dimension 6. In the spectral function we take the p meson in the narrow width approximation and the 
^ '. NN stated 

A complete or "full" in-medium sum rule is actually the corresponding vacuum sum rule, perturbed by small, density 
dependent terms. Such a "full" sum rule, when evaluated, is guaranteed to be stable, more or less, irrespective of 
these small contributions, simply because the vacuum sum rule is so. To get rid of this insensitivity, we subtract out 
the vacuum sum rule from the "full" one, retaining on both sides only the terms to first order in the density. 

Though well-known, these sum rules are derived here to bring out some technical points. We start with a broader 
framework, namely the ensemble average of the correlation function of currents at finite temperature and chemical 
potential. The special cases of sum rules at finite temperature and at finite density may then be recovered from the 
general formulation by setting respectively the chemical potential and the temperature to zero. We derive the mixing 
of two operators of dimension 4, namely the quark and the gluon parts of the energy momentum tensor, a result 
known in the context of deep inelastic scattering. We also present a field theoretic derivation for the NN contribution 
to the spectral function. 



*Electronic address: mallik@tnp.saha.ernet.in 
^Electronic address: nyffeler@cpt.univ-mrs.fr 

'Actually there could be further density independent contributions in the continuum, but they are eliminated in the subtracted 
sum rules we shall be considering. 
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On evaluation we find that the sum rules as saturated above do not yield stable results, as we have already reported 
recently |J . This instability points to the importance of other contributions we did not include, namely, the Lorentz 
nonscalar operators at dimension 6 and the Landau singularities besides the one from the NN state. The neglect of 
the (large) in-medium width of the p-meson also possibly adds to this instability . 

In Sec. H we introduce the kinematic decomposition of the two point function of vector currents. In Sec 



[II a 



simple derivation of the operator product expansion is presented in coordinate space. Then in Sec. IV we construct 
the spectral representation from p and NN intermediate states at finite density. The sum rules are written and 
evaluated in Sec. |v|. We conclude in Sec. VI. 



II. KINEMATICS 

Consider the ensemble average of the time ordered (T) product of two currents 

T$(q) = i J d'xe^ (T{V«(x)V»(0)}) . (2.1) 
Here V®(x) is the vector current (in the p-meson channel) in QCD, 

Vftx) = q(x) 7 ^q(x), (2.2) 

q(x) being the field of the u and d quark doublet and r a the Pauli matrices. The ensemble average of an operator O 
is denoted by (O), 

(O) = Tr \e-^ H -^0] / Tr Lh^-mAO] , (2.3) 

where H is the QCD Hamiltonian, f3 is the inverse of the temperature T, and Tr denotes the trace over any complete 
set of states. Keeping our application in mind, we have introduced the nucleon chemical potential \i with the 
corresponding number operator N. 

For explicit calculations one generally chooses the rest frame of the medium, where the temperature is defined. This 
breakdown of Lorentz invariance may be restored if we let the medium have an arbitrary four- velocity []l2"|| . [In 
the matter rest frame u M = (1,0,0,0).] The time and space components of g M are then raised to the Lorentz scalars, 
w = u ■ q and q — y 'w 2 — q 2 . 

There are two independent, conserved kinematic covariants |l3| ], in which the two point function from Eq. (2.1) 



may be decomposed. In choosing their forms we must note that the dynamical singularities extend up to q 2 — 0. We 
thus have to ensure that the kinematic covariants do not have any singularity in this region, particularly at q 2 = 0. 
The covariants P^ v and Q^ v defined by 

p q^_^_~ . n _q*_~ . , 9 .s 

q A q z q z 

where — — wq^/q 2 , are indeed so, as can be seen by inspecting their components for space and time indices. 
The invariant decomposition of is then given by 

T$(q) = 5 ab (Q flv T l + P^T t ), (2.5) 

where the invariant amplitudes Ti t t are functions of q 2 and w. 

The kinematic covariants are still not regular: they have a dependence on the direction of the three-vector q even 
as \q\ — > 0. Thus the spatial components of Tff in this limit become 

T«\q) = S ab [S l3 T t - mnj(T t - qlm, (2.6) 

where ni is the ith component of the unit vector along q. To eliminate this direction dependence, we must require 
that § 

Tt(qo, \q\ = 0) = ggT,(9b, l?1 = 0). (2.7) 
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III. OPERATOR PRODUCT EXPANSION 



Here we outline the derivation of the short distance expansion of the product of two current operators in Eq. (2.1). As 
already stated, the availability of the velocity four-vector of the medium allows one to construct new, independent 
Lorentz scalar operators. Thus up to dimension 4 we have, in addition to the old ones 1, qq — uu + dd, G 2 = 
{a s j-K)G a illv G lll ' a , appearing in the vacuum expectation value, two new ones Q q = u^O^u" and Q 9 = u^Q 3 u v for 
the ensemble average. Here G c v , c = 1, 8 are the gluon field strengths and a s = <?f/47T, g s being the QCD coupling 
constant. are the quark and the gluon parts of the traceless energy momentum tensor 

Tfi 

®liv — ~G € ^ x G v Xc + -g f _ l „G c af3 G al3c , (3.1) 



where m is the quark mass in the limit of SU(2) symmetry. 

We now obtain the (singular) coefficients of the operators in coordinate space jl4|-|l7j . In the free field theory these 
are obtained from the Wick decomposition of the operator product. If we treat the gauge fields as external, the same 
form of decomposition also holds for the interacting theory 

T {V°(x)V„ a (0)} = tr[^S(x, 0)7,5(0, x)r a r b /4] 

-iq(0)j„S(0, x)(T b T a /4)q(x) - iq(x) 7fl S(x, 0) 7 „(tV 6 /4)g(0) 

+a regular term. (3-2) 

The trace in the first term is over all the indices, namely spin, flavor and color of the quark field. The quark propagator 
S(x,0) is now in the presence of the (matrix valued) gauge potential Au(x), 

{-i-f[d M - igA M (x)] + m) S(x, 0) = 5\x). (3.3) 

An important technical point here is that quantities can be put in a gauge covariant form from the beginning, if 
we use the Fock-Schwinger gauge, in which, x M Ap(x) — 0. In this gauge 

A^x) = (l/2)x Q G QM (0) + --- , 
q(x) =q(0) +x»D tl q(Q) + ■■■ . (3.4) 



Equation (3.3) may be solved for the quark propagator in a series of increasing number of gluon field strengths. 

The unit and the quark operators 1, qq, and Q q reside in the first and the next two terms in Eq. ( [3.2] ), 
respectively. To calculate their coefficients, it suffices to use the free quark propagator 

„ , If 2af irh \ 

S (x, 0) = -j - 2 * - , (3.5) 



47T 2 \ (x 2 - ie) 2 



where we retain only the leading singular term in rh. The quark operators in Eq. (3.2) can then be projected on to 
the scalar operators by using the formulae 

(VaQb) = l^ k (S BA (qq) + (+) BA {ftq)), (3.6) 
(q A i{D,q) B ) = ±5* \^5 BA u,m(qiq) + (j,) BA (jt(qq) - i<0«>) + ^(^^(9')} . (3.7) 

The flavor and spinor indices of the quark field have been denoted by (j, k) and (A, B), respectively. A sum over color 
indices on both sides is understood. These projection formulae are easily obtained by expanding the quark bilinear 
matrix in terms of the complete set of Dirac matrices, noting parity conservation, and using the free equation of 
motion for the quark field. (Chirality forbids (qtf.q) to appear in the vector-ve ctor correlation function.) 



The two-gluon operators G 2 and <d 9 arise only from the first term in Eq. (3.2). To get their coefficients we have 
to work out the quark propagator up to two gluons. The general fourth rank tensor with two gluon fields may be 
projected on to the scalar operators by 
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(t? c G aP G Xa ) = -(g a xgp. - g a ,gpx) - A {G"G^ C ) + (6«) 



(3.8) 



tr c indicating trace over the color matrices. It turns out that if we insert this expression into the two gluon terms 
present in S, they do not contribute to G 2 . Thus the calculation of its coefficient reduces to inserting S up to the 
one-gluon term 



S(x,0) = S (x,0) 



16tt 2 (x 2 — ie) ' 



(3.9) 



in the first term of Eq. (3.2) and then using Eq. (3.8). However the two-gluon terms in S do contribute to Q 9 . But it 
is not necessary to calculate this coefficient. As we show below, the renormalization group improvement automatically 
brings it in. 

Having indicated the procedure to calculate the coefficients in coordinate space, we take the Fourier transform and 
write the result in momentum space (Q 2 = —q 2 ), 



i I d 4 xe lqx T{V^{x)V*(0)} = S ab 



\8tt 2 



1 



-L H Q 2 /^ 2 )l + ± 



1 



1 



-mqq 



24 



:G 2 



2? 
Q 2 



-mqq 



0« 



—g 2 + -e« 

24 3 



(3.10) 



Here (i is the renormalization scale, which we take to be 1 GeV as the natural scale to normalize the operators. But 
this scale is not convenient for the coefficients. These, calculated above only to the lowest order in a s , may contain 
powers of \n(Q 2 / /i 2 ) in higher orders. To get rid of these large logarithms, |j we must shift the renormalization point 
for the coefficients from /i 2 to Q 2 . 

Now the left hand side of Eq. (3.10) is independent of the renormalization scale. So are all the operators 1, mqq, G 2 
and hence their coefficients. But it is not so for 8 9 and 9 .[] Indeed, it is well-known in the context of deep inelastic 
scattering that 6 9 and 6 9 mix under a change of scale . Let us denote the contribution of these operators generally 
as 



Gq |p <d v |p 



~G g \n O 9 |p= ^ Ci \n 0' | M , 



(3.11) 



where only the dependence on the scale fi is shown explicitly and the index i runs over q and g. The coefficients Ci 
satisfy a coupled set of renormalization group equations having the solution 



Gi Cj 



- r o S(t')di' 



(3.12) 



J1 



where t = 4 ln(Q 2 / /i 2 ) and 7 is the 2 (g> 2 anomalous dimension matrix. In the basis we are using it is 



7 = 



(4tt) 



3; / 



2 l M _ 



i 



3 U f 



(3.13) 



where «/(= 2) is the number of flavors. The exponentiation in Eq. ( 3.12) can be easily worked out by a spectral 
decomposition of the matrix 7 = J^. Ai-Pj, where Ai and Pi are the eigenvalues and the corresponding projection 
operators p9|. The eigenvalues are 



Ai = 0, A 2 



4/16 



(4tt) 2 3 



; 7 + n f 



(3.14) 



2 The logarithm in the coefficient of the unit operator in Eq. (3.10) is of different origin. See the first footnote in Sec. 4.6 of 
Ref. §. 

3 If we had calculated the coefficient of G 9 , we would have found it to be proportional to a s ln(Q 2 //i 2 ). It signals mixing with 
other operators. 
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and the projection operators can be built out of the (right) eigenvectors and left eigenvectors. We thus finally get 



(¥+»/; 



„ 16^, 



16 

+ A(Q 2 )(C 9 -C 9 ) | Q ( 



where = 9 9 



s and the anomalous dimension of the second operator resides in X(Q 2 ), 

-d/2b 

HQ 2 ) 



a s (ii 2 ) 



a s (Q 2 



d = 



16 

t 



71/ 



b=ll- -Uf. 



(3.15) 



(3.16) 



We will use a s (/i 2 = 1 GeV 2 ) = 0.4 in our numerical evaluations below. 

Once the renormalization group improvement is made, we can replace the coefficient functions by their lowest order 
values. These may be read off from Eq. (3.1C) for the amplitude Tj, for example, as Cf \q= 2/(3Q 4 ), C g \q= 0. Thus 
the renormalization group improved result for XJ reads (with nf = 2) 



24 



HQ 2 



3 09 



e 9 



(3.17) 



Notice that if X(Q 2 ) is set equal to unity, it collapses to the previous result in Eq. ( 3.10| ). The amplitude T t is also 
given by the same expression, except for an overall factor of — Q 2 and a factor of (1 — 2q 2 /Q 2 ) multiplying the term 
with <d's. 

Let us summarize the effect of mixing of g and 9 under the renormalization group. It introduces the multiplicative 
factor, X(Q 2 ), which is nonperturbative in the sense that it is a sum over all the leading logarithms in the perturbation 
expansion. It is true that the coefficient of 9 is smaller than 9 by a factor of g 2 . We have left out this smaller 
coefficient, as did the earlier authors, but retained the nonperturbative factor. Numerically the factor is significant 
for Q 2 away from p 2 = 1 GeV 2 . 

At dimension 6 we include only the Lorentz scalar operators in vacuum, namely the four-quark operators |EJ. In 
the approximation of ground state saturation, its contribution reduces to the familiar expression involving (qq) . 



IV. SPECTRAL REPRESENTATION 

QCD sum rules make contact with the physical states through the spectral representation of the two point function. 
For the vacuum amplitudes, these are the Kallen-Lehmann representations in q 2 , the only variable available there. 
At finite temperature and chemical potential, it is convenient to use the Landau representation, which is a spectral 
decomposition in qo a t fixed \q\ pG]. For the invariant amplitudes they are 



1 f + °° 
Ti,t(qo, \q\) = - / dq 
n J-oo 



ImT i,t(gd ; Iffl) 
q' -Qo- ie 



tanh(/3^/2), 



(4.1) 



up to subtraction terms. The imaginary parts from different intermediate states can be written as phase space integrals 
over the matrix elements of currents. But we shall not write these expressions directly. Instead, we first calculate the 
full amplitudes with these intermediate states in the real time formulation of field theory at finite temperature and 
density |^l| and then read off their imaginary parts. Since we have to calculate the amplitudes at the tree level, only 
the 11-component of the propagators will be needed. 

We first consider the p-meson pole, see Fig. 0(a), which would dominate the absorptive part of the vacuum correlation 
function. 




(a) (b) 

FIG. 1. Contributions to the spectral function from (a) the p meson and (b) the AW state. 
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In the vacuum, its coupling to the current V p (x) is given by 

(0\V«(0)\p b )=6 ab F p m p e p , (4.2) 

where e M and m p are the polarization vector and the mass of the p. Experimentally F p = 153.5 MeV. This matrix 
element suggests that, as far as the p-meson contribution is concerned, we may write the operator relation 

V«{x)=m p F p pl{x). (4.3) 

In the medium such a relation is generally modified in two ways. First, the parameters m p and F p must be replaced 
by their effective values in the medium, to be denoted below by a star. Secondly, the space and the time components 
have different constants of proportionality. However, the latter modification tak es place in higher orders in density 
and temperature p2fl . As we work to the lowest order, the operator relation ( [1.3] ) changes in the medium simply to 

VZ{x) = m*F* pP l{x). (4.4) 

Then the p-meson contribution to the correlation function is given essentially by its propagator in the medium. As 
before, we get 



I ! 1 1 '/ } ( C/.; ) . if 1 ] \ / 1 ^ _,.■■>.<. ,|,7l-> , 

lmT t (q , \q\ 2 ) J { m* p 



\ 2 nFfd (q 2 - (\q\ 2 + mf)) coth(/3g /2). (4.5) 



In the medium there are other equally important contributions. These are from the 2ir and, if the medium is nuclear, 
also the NN intermediate state, see Fig. [l](b). As we restrict in this work to sum rules at finite nuclear density but 
at zero temperature, the 2tt contribution would be small and close to its vacuum value. The NN contribution in 
vacuum would also be small, beginning with a cut at the threshold q$ = Am 2 N + \q\ 2 . However, in the nuclear medium 
the current also interacts with the real nucleons to give rise to a short cut around the origin, — \q\ < q < +\q\. Below 
we evaluate this contribution. 

The contribution of the nucleon isodoublet, spinor field N\(x), i = 1,2, A = 1, 4, to the vector current is 



V£(x) = N(x) ltl T -N{x). (4.6) 



The NN contribution to the two point function is then 



Vfa) = 2 6a J ^xe^tri^SixM-x)}, (4.7) 

where the trace is over the gamma matrices and iS(x) is the 11-component of the nucleon propagator in the medium 
(T {N\(x)N' B (y)}) = iS»S AB (x - y), with 

iS{x-y)= [ ^L e - ik < x - y \^ + m N )\- 2 V— --^5{k 2 -m 2 N ){n-9{k )+n + e{~k )}\, (4.8) 

J s nijj + if | 

the distribution functions for the nucleon and the anti-nucleon being given by ri*{u) = (e^TP) + w ith 



UJ 



'N- 



Our aim is to integrate out the time component of the loop momentum, when the imaginary part of T p b (q) can be 
read off from the resulting expression. However, the above expression for the propagator is not convenient for this 
purpose because of the presence of 8(±ka). So we carry out the fco integration in the propagator itself, getting 

/r1 3 k 
[eix^i^ + m N )(l n-)e- ik -* - + m N )n+e^} 

+ 6(-x°) {-(ft + m, N )n-e- lk - x + + m N ){l - n+)e lk - x }] . (4.9) 



Here the time component of fc M is understood to be given by u>. We can now integrate over x^ in Eq. (4.7) to get 
(5-functions in three-momentum and energy denominators in the time components. After some algebra we get for the 
imaginary part 



G 



where 



with 



ImT^(g) = --S ab coth(f3q Q /2){I^(q)+I^(-q)}, 



d 3 k A M „ 
(27r) 3 4cjiw 2 



[(1 - n 1 - n%)5(q - u>i - lo 2 ) - (n 2 - n 1 )S(q - uji + w 2 )], 



(4.10) 



(4.11) 



A M „ = tr{7 M (^ + m N )~i v (ty - fa + m N )}, 

= AftkftK ~ (k^q v + Kq^) - g^{k ■ (k - q) - m 2 N }}. 



(4.12) 



Here the subscripts on n T serve to indicate that their arguments are W\ = y k 2 + m 2 N and lu 2 = y (k — q) 2 + m 2 Nl 



respectively. The first and the second terms in Eq. (4.11) are nonvanishing in the timelike (q 2 > 4m^) and spacelike 
(q 2 < 0) regions, respectively. Working out the angular integration, it becomes 



V(<?) 



1 



duii 



i\q\ J (2tt) 2 ' A ^ 1 ~ 711 ~ n2 "> 6 ( q2 ~ 4m ^ 



(4.13) 



where the integration limits are given by uj± — [q a ± \q\v(q 2 )]/2, with v(q 2 ) = yjl — 4m^ / q 2 . 

To project the tensor on the imaginary parts of the invariant amplitudes it is convenient to form the scalars Ifi 
and u^I^u 1 ' and relate them to the invariant amplitudes. Changing the variable oj\ to x given by uii = ^(qo + 
we finally get for the contribution from the NN state 



ImT,(go,M 2 ) 
lmT t {q Q ,\qf) 



coth(/3g /2) 



v{i-v 2 ) ( 1 
24tt 2 



for q 2 > Am 2 N , 



(4.14) 



with 



If 
It 



32tt 2 



dx 



2{x 2 - 1) 



(2 - v 2 + x 2 ) ) {n ^ qo + l^/ 2 n ^ qo ~ ^ x)/2 + 



(4.15) 



and 



ImT ; (g ,kT) 
JmT t (q , \q\ 2 ) 



~coth(0q o /2)( ^_ 



for q 2 < 0, 



(4.16) 



with 







a: 


I 32** / 



dx 



2(x 2 - 1) 
-q 2 {2 - v 2 +x 2 ) 



2(n „2\ji\ ) {n((-q +\q\x)/2- fj,) -n((q + \q\x)/2- fi)}. 



(4.17) 



The superscripts (±) on // ;t denote timelike and spacelike q^, respectively. The distribution function n appearing in 
Eqs. fl4T5| ) and $~T% is defined by n(e) = (e Pe + l)" 1 . 

Still higher continuum contributions, at least their leading parts, are density independent, as is shown by the 
coefficient of the unit operator. They will drop out in our subtraction process. 



V. SUM RULES 



We now go to the spacelike region in q^, (Qq = —q 2 > 0) at fixed |g| and take the Borel transform of both the 
spectral representation and the operator product expansion with respect to Qq. The sum rules are obtained by 
equating these transforms at sufficiently high M 2 . For TJ it is 
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F *2 e -m;>/M* + 1 / dq ^ /M^3 {q 2 + |gf) 

+ e \e/M> ( r dgge-^I+fo, m + f m dq l e -^l- (qo , m ) 

\Ji m l + \q\ 2 JO J 

M 2 (O) (0>) 

8ir 2 M 2 2A/ 4 ' [ ' 

where (O) and {O') represent the contributions of all dimension 4 operators and of dimension 6, scalar operators in 
vacuum, respectively, 

(O) = \m{qq) + ^ + ^ {(6) + A(M 2 ) (!±Q* 9^} , (O') = ^(m) 2 - (5.2) 

The amplitude T t satisfies a similar sum rule. 

Considerable simplification results if we take \q\ — > 0. This limit may be taken inside the integrals in Eq. ( f5.l| ), 
except for the last one, where the integrand diverges, while the range of integration shrinks to zero. A careful 
evaluation gives a result resembling a pole at q 2 = 0.^ Then the above sum rule for T; becomes 

1 r°° r 

F *2 e -m* p /M* + __ I dq 2 Vo{3 _ v 2 q) \ n{go/2 _^ )+ e ^ 2 o/M 2 {1 _ n{qo/2 _ /i} _ n{go/2 + M)} 



M 2 (O) (O') 



(5.3) 



8tt 2 M 2 2M 4 ' 
while the sum rule for T t in this limit may similarly be found to be 

1 r°° 

m *2 F *2 e - m ?/M* + dq 2 q 2 v (3 - ^e"^ 2 {1 - n (q /2 - M ) - n(q Q /2 + M )} 

-5i-OT + y. <M> 

where «o = (1 — 4to^/(/q) 1 / 2 . Note that the contribution of the short cut is nonvanishing in this limit only in Tj. 

As they stand, the above sum rules cannot be well saturated, since the 2ir and the continuum contributions on the 
spectral side are missing. However, if we subtract the corresponding vacuum sum rules, the continuum contribution 
will practically drop out. We now go to the limit of zero temperature, when the 2tt contribution reduces to zero. 
With [i > 0, the nucleon and the antinucleon distribution functions reduce in this limit to a ^-function and zero, 
respectively. We thus finally arrive at the following sum rules at finite nuclear density: 



Ffe'-f'^ +^f] ds(l - e-^^l-H/^l + 2m%/ S ) = g> - g> , (5.5) 
m ?F; 2 e- m ?' M2 - r ds a e- a ' M2 yjl 4m 2 N /s(l + 2m%/s) = -(O) + -gj, (5.6) 



where the bar over a quantity denotes subtraction of its vacuum value, e.g., (O) = (O) — (0|O|0). 

The spectral integrals in the sum rules can be expanded in powers of the Fermi momentum pf = \[W 
the nucleon number density W 



m 2 N , or 



' /! V(PF-b1) = M- (5-7) 



(2tt) 3 ^ nu in 2 ' 

The operator matrix elements can also be expanded in powers of the nucleon number density. With the normalization 
of the one nucleon state (p, a\p' , a') = (27r) 3 2po<5a,a"5 3 (p — p*), a, a' denoting the spin and isospin degrees of freedom, 
we have the expansion for any operator R up to first order as 



4 The situation is similar to that at finite temperature Q. See also Ref. 
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FIG. 2. Coefficients of n/n a in the expansion of m* and F* from Eq. (5.17) as function of the Borel parameter M 



A numerical evaluation of a and b is shown in Fig. |^. It will be noted that there is no sign of constancy of a and b in 
any region of AT 2 ; rather, they increase rapidly and monotonically throughout, in particular a. We conclude that the 
sum rules as saturated above are inconsistent and cannot give any reliable information about the density dependence 
of m p and F p . 

VI. DISCUSSION 

We have studied QCD sum rules in the nuclear medium in the p-meson channel using a saturation scheme considered 
earlier by a number of authors The sum rules we actually evaluate consist of terms, all proportional to the 

nuclear density. We find that the results are not stable in any range of variation of the Borel parameter. We 
conclude that the sum rules with the present saturation scheme become inconsistent, in apparent contradiction with 
the evaluation by these authors. 

The resolution of this contradiction lies in the fact that the "full" sum rules considered by these authors are 
dominated by vacuum contributions, the density dependent terms serving as small perturbations. As stated already 
in the Introduction, their stability follows essentially from that of the vacuum sum rules. The variation of some free 
parameters such as so, the beginning of the continuum, in these works would further contribute to this apparent 
stability. 

It would appear from earlier works that QCD sum rules can accommodate different saturation schemes, each one 



giving rise to a new set of results |11|. By working with a more sensitive version of the sum rules, we show that the 
situation may not actually be so: a simple saturation scheme, such as here, is rejected and a more realistic one is 
called for. 

Let us finally discuss the neglected contributions, whose inclusion should restore the stability of these sum rules. At 
dimension 6, five more operators would have contributed. These contributions can be evaluated using the data on deep 
inelastic scattering ^6|. To get an idea of the convergence of the series of operators, we note that the contributions 
of a family of similar operators, differing only in the number of Lorentz indices, form a series in powers of (mjsr/M) 
with decreasing coefficients, given by the corresponding moments of the structure functions. 

In the spectral function, it is not just the NN state, which can give rise to the Landau cut in the low q 2 region. 
Any other state NN*, where A* is a resonance communicating with the pN channel, can contribute a cut for 
q 2 < (tojv* — tin) 2 - The two resonances A(1232) and N(1440) would presumably contribute significantly. To include 
these contributions we need, however, know the current- A A* couplings. Also the increased width of the p meson in 
the medium is another source of altered contributions to the sum rules. We note that a study of such "subtracted" 
sum rules for p mesons at finite temperature showed no such instabilities as mentioned above |27| . 
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